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Abstract 

Probabilistic estimates on linear combinations of eigenvalues of the one dimensional Anderson 
model are derived. So far only estimates on the density of eigenvalues and of pairs were found 
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linear Schrodinger Equation, where linear combinations of eigenvalues are the denominators and 
evaluation of their smallness is crucial. 
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I. INTRODUCTION 



The study of the statistical properties of eigenvalues and their distribution for the An- 
derson model has attracted a lot of attention in the Mathematical Physics and Physics 
communities. The one dimensional Anderson model on a lattice can be defined as an eigen- 
value problem 

E^ E (x) = ^ E (x-l)+^ E (x + l)+e x ^ E (x) xeZ (1) 

where ip E {%) is the eigenfunction corresponding to the eigenvalue E. The sites of the lattice 
are denoted by x while e x are independed random variables which are unifomly distributed 
in the interval [—A, A]. It exhibits compelete localization and does not have any exceptional 
states. Important progress in the understanding of the distribution of eigenvalues has been 
made at the rigorous level with the fundamental works of Wegner [1] (also see e.g. Aizenman- 
Molchanov (AM) fl) and Minami [j], U], which describes the main aspects of the distribution 
of the distances between the eigenvalues. Wegner's estimate shows that the typical distance 
is of order l/V, where V is the volume of the cutoff Hamiltonian. A lot of numerical 
and other arguments, led to the conjecture that the distribution of the distances between 
eigenvalues (EV) comes as from a Poissonian distribution, and the probability to have two 
EVs in an interval of size I is proportional to I 2 . This last statement was finally proved 
by Minami Q|. Very recently, we showed that the distribution of the energy level distances, 
and higher order combinations are crucial to understanding of the large time behavior of 
;he Nonlinear Schrodinger Equation (NLSE) with the Anderson potential in its linear part 
l]]. This can be understood when one studies the appropriate time dependent perturbative 
construction of the solution, around the linear part. Then, the leading potentially large 
terms consist of denominators which are linear combinations (with integer coefficients) of 
the various EVs of the linear or renormalized linear problem. Thus we are naturally led to 
consider the probabilistic estimates on the smallness of quantities like 

R 

f = E c ^ ( 2 ) 

k=l 

where the sum over k is finite. Here are integers and E^ are the eigenvalues of the Anderson 
model (pQ). The indexing of the EV we use is such that to Ei corresponds the eigenfunction 
which is localized around the lattice site i 



■ in dul 



Very recently progress on this problem has been achieved in the works j|| and Q|. Our aim 
in this work is to obtain the needed estimates for the perturbative construction of large time 
solutions for the NLSE with Anderson potentials jjL(j| , as developed in [3, [ill] - Our strategy 
is to begin with the AM approach [2}], and estimate the expectations of fractional powers of 
the inverse of the linear combinations of EVs, denoted by / (Sec. III). It turns out that the 
critical step in this problem is to estimate from above the probability of avoided crossings, 
which is equivalent to controlling the probability of states, with two humps that are far 
away. Another useful tool is the demonstration of a lower bound on the eigenfunctions, 
which is uniform in the potential, with large probability (Sec. II). A modified version of / 
is considered in Sec. IV. 



II. THE ASYMPTOTIC EXPONENTIAL DECAY OF THE EIGENSTATES 

For each energy E the inverse localization length, 7 (E) , is continuous as a result of the 



Thouless formula |12j] and, 7 (E) > (see more details after (JTTJ) ) . Since the spectrum of a 
Hamiltonian H on a lattice is compact, it follows that, given f), we can cover it with a finite 
collection of open intervals, such that each interval is around A&, k = 1, 2...K < 00, and for 
energies in that interval (1 — 77) 7 (A*.) < 7 (E) < (1 + 77) 7 (A^) , where f] is arbitrarily small. 

Now we focus on one of such intervals located around the energy A 6 J , and fixed 
e <C 2f/7 (A ) • We denote by s n j the set of all the potentials for which the energy Ej e Iq, 
and such that, 

| e (7(%)-H0M^ I > 1 je z (3) 

the inequality holds for all |x| > n (fj,e) , and x is the distance from the localization center 
of the wave function. It is a consequence of the Furstenberg theorem [13j| and of the analysis 
of that n is finite for a.e. potential (for an alternative approach see |l5|). 

The measure of the potentials belonging to s n ,j is denoted by S n j. Clearly, S n ,j > 0, and 
we only consider the case where the measure of the set of potentials with Ej belonging to 
the interval around Ao is nonzero. S n j is a monotonic increasing sequence of n: this follows 
since by construction, if S n j is such that condition ([3]) is satisfied for \x\ > rij it is also 
satisfied for \x\ > rij + I, I > 0. Furthermore, one can verify that each S n j is measurable. 

Lemma 1 The limit Hindoo S n j = fi' Q > 0, exists with some positive fi' , where fi' is the 
measure of all the potentials, such that Ej is in the interval around Aq. 



3 roof. The lemma follows, since by Furstenberg theorem [13| and Delyon-Souillard-Levy 
for almost all energies the eigenfunctions decay exponentially with the rate 7 (E) : 



In 



lim 

x— >oo 



1 1 2 1 1 2 



2x 



7 (Ej) 



Therefore 



1 /2 

^ (x)\ 2 +\^ Ej (x + 1)| 2 ) x P(x)e-^W 



(4) 



(5) 



with P (2) sub-exponentially bounded, above and below (decaying or increasing), and x is 



the distance from the localization center 116 
Hence 



121 



lim \e^ )+£) ^P(x)e-^ E ^\=oc 



(6) 



for all £ > 0, and therefore for each eigenfunction condition j3]) is satisfied for some finite rij. 
Hence all realizations of the potential satisfy condition j3]) for any E in the interval around 
A : 

fio = fi (all potentials) = fi [Js nd < ^/x (s n>j ) , (7) 



11 ,3 



since {n, j} is countable. ■ 

The lemma holds for each of the open sets around Since the number of these sets 
is finite, then for any fixed rij the condition (J3]) is satisfied for a set of realizations with a 
measure of 1 — 8j, where Sj decreases with rij. 

This leads to the following proposition. 

Proposition 2 (Lower bound on Eigenfunctions) Given e,5 > 0, with probability 1 — 8, 
3 {n* (e,6)} C N, such that 

1 /2 

\U) Ei {n)\ 2 + \tp Ej (n + l)f) x P(n)e-^ )W > e~ (7( ^' )+£)W (8) 

/or all n > n* (e, 5), and where n is distance from the localization center. The proof of the 
upper bound is similar. 

R 

Corollary 3 For any f = '^^CkE ik , there is an x such that for j > \x\ , 



k=l 

Pr 



df + df 



<C 3 )<6(j) 



(9) 



dEj dej + i 

where Cj > e^ 2 ^ +£ ^ and 5 (j) is monotonically decreasing with j. The decay rate 7 corre 
sponds to an energy E ik in the sum for f , namely 7 = ^y(E ik ). 



Proof. Using Feynman-Hellman theorem 

df + df R 



dej de 



(10) 



5>0^ m 2 + \A k (j + i)i 2 ) 

k=l 

First we order the eigenfunctions in the RHS of f TTOT ) by their localization centers, 
{%i < %2 < ■ ■ ■ < in} , where R is the number of eigenfunctions in the sum. If the eigen- 
functions are tightly packed in the sense, that their localization centers are found in a 
box of size n* (e,5) then we can take j such that \j — ir\ > n* (s,8) and Proposition [2] 
would apply for all the eigenfunctions in the sum. We assume for simplicity j > 0. Here 
n* (e, 5) = sup k { im jfl ^ (e) . For j large enough the asymptotic behavior of this sum is 
dictated by one term with the smallest inverse localization length, 7. In the worst case, 
there are two states which are localized at a distance of n* apart on the lattice and have 
very close inverse localization lengths 7, (j±r]), with r\ small. By continuity of 7 [E] the 
respective energies are typically close. For these states we have to go for a distance of 

j > lmmn*/r], (11) 

where 7 min = inf^ (E) . This results from the requirement that e~ lrai ^ ^> e~^ iTaXaJrrt ^~ n > , 
which is satisfied for e^ min+ '^ n * ^ 1 . We will disregard realizations of the potentials for 
which many 7's are closer than 77, however we notice that the measure of such realizations 
goes to zero with rj. In fact the smallest energy spacing results from avoided crossings which 
are exponentially small in n*. In this process we practically eliminated realizations where 
arbitrarily small energy differences are found. (See next section for detailed discussion of 
this issue, a quantitative estimate is (|T9l) ). For a general configuration of the localization 
centers in the sum we will separate the localization centers to clusters such that the distance 
between two adjacent clusters is at least n*7 min /^. Following from above, there are intervals 
on the lattice between the clusters, such that the contribution of each cluster to the sum is 
governed by one exponential. Therefore we can always find j's between the clusters such 
that only one cluster and actually one state dominates, moreover at this point the magnitude 
of the sum satisfies the inequality 

Cj > e" 7 — 7 --™*^, (12) 
where 7 max = sup E 7 (E) . In this way we obtain the following proposition. 
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Proposition 4 Given R energies E^, Ei 2 , Ei R we can always decompose this set into a 
set of disjoint clusters such that the largest n* is cR/^ m i n , where with probability 1 — e, for 
j > n* (e, R) 



k=l 



> ce 



-2(7+e)b'l 



with Cfc that are integers and c is a constant. That is for some j , some 

2 

all other ipEi (j) 



dominates 



Using the above proposition we can find an x, such that for j > x, the sum (jTOJ) will be 
dominated by one term with the smallest 7 (E) , such that 

Of 



Pr 



<CA <5(j) 



(13) 



and Cj depends only on the relative distance between j and the localization center of the 
dominant eigenfunction, and therefore does not depend on the locations of the other local- 
ization centers of the states in /. Here 5 (j) is the measure of all the realizations that were 
excluded so that Proposition [2] holds. Due to the Proposition [2] both Cj and 5 (j) , confined 
with fpTTT ) are decreasing functions of j. Therefore Cj and 5 (j) are bounded by C x and 5 (x) 



III. THE INTEGRAL OVER \f\' s 

Lemma 5 For each j, exists a finite interval Asj > for which the function Jj- + d ^! +l 
does not change sign with probability 1 — Si(j) 



Proof. The function 4^- + a 9 ^ is given by 



j+i 



df_ , df 

dsj 



+ 



de j+1 



R 
k=l 



U) + ^ t U + 1) 



(14) 



Due to Corollary [3] there is a j > jo for which the sum is dominated by only one term with 
the smallest 7 (E) . We will denote the state corresponding to this energy by its eigenvalue 
Ej. The sign of -f^ + can change only when there is an avoided crossing between the 



Ej state and one of the other states Ek, in the sum defining /, since then the dominating 
terms will interchange due to the fact that 7 (E) is a continuous function of E. Consider a 
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cluster of size L, which is large enough to cover all the eigenstates, such that outside of this 
cluster the relevant eigenstates are smaller than e~ 7minn *. Using this cluster, the Minami 
estimate gives 

PrfTrP£ ) (/)>2)<( 7 r||p|| oo /L) 2 , (15) 



where WpW^ is the supremum of the density of states, / is some energy interval while Pjj\l) 
is the spectral projection on that interval. We can reformulate this estimate in such a manner 
that the probability to find any i and k with energies in this interval satisfies 

Pr(| J E i - J E fc |</)<(7r||p|| oo /L) 2 . (16) 

Dividing the spectrum into intervals of length / = b\L~ v with b\ > and v > 2, we find 
that the probability to find any pair i, k in such an interval satisfies 

Pr {\Ei -E k \< hL-») < ^Att 2 \\p\\l L~^~ 2 \ (17) 

Therefore the probability of an energy interchange can be made arbitrarily small. After the 
interchange of the energies the magnitudes of the corresponding 7's will interchange and one 
(different) state will dominate at j. This assumes that d^/dE for this state is not too small. 
If it vanishes the following evaluation is required. In this situation, the change in 7 for a small 
change £3 in E, we assume that there is a finite number of points, Nb, where the derivative 



of j(E) is zero, w 



lich follows from the analyticity of 7 (E) , that in turn follows from the 



Thouless formula l^] and the analyticity of the density of states [if], Q, 0|- We take a 
small interval, 5 2 , around those points for which, d^/dE < 5 2 . We will exclude realizations 
with energies which are found within this interval, by using the Minami estimate, 

Pr {\E, -E k \< I(S 2 )) < N b (tt HpII^ I(5 2 )L) 2 < ^(j). (18) 

where 1(62) is the interval of energies for which , d^/dE < 5 2 - This interval, /(<5 2 ), is a 
monotonically decreasing function of 5 2 . Energies which are found outside of this interval 
will have a derivative d^/dE > 5 2 and therefore using (fPTj) 

Pr (MEt - 7(^)1 < tfuWr") < tt^A HpHL L- (2 "~ 2) < Uj)- (19) 

where 5 2 (j) can be made arbitrarily small. Here we choose L to be much larger than the 
cluster in question. To ensure that this probability is small we have to choose v > 2. 
Note that the number of intervals grows as L v . The realizations which are excluded are of 
probability which is bounded by Si(j) = 6i(j) + ■ 



Remark 6 The bound on the probability of the previous corollary should be modified with 
5j replaced by 5(j) + 8i(j)- 

R 

Theorem 7 For f = ^^CkEi k ,the following mean is bounded from above 



k=l 



1 

17iY> 





A 


h 


= 1 



d\x (ei 



\f\ l 



(20) 



where (.) s denotes a mean over realizations of a measure 1 — 5 and D$ is independent of 
|A| .We analyze the behavior of the integral 

A 



'/ 



d\i (ei 



1 



i=i 



\fm 



(21) 



Proof. We denote, e = {ei, • • • £|A|} , and assume, that is uniform in the interval 

[—A, A]. We first change variables (for a j chosen to optimize the bound on (122]) ) to 

i 

with the Jacobian equal to 1. Then we change the variables to (/, S\ 
that produces the Jacobian 



(ej + e j+1 ) 

[£i — £j+i) 



: |A| 



)■ 



d(f,e 1 ,...,e lM ) 
d (ei, . . . ,£| A |) 



df df 
dsf dei 



df 



|A| 



df 



del 



(22) 



The integral using the new variables is 



(2AfJ d£j J-A 



MA) 
dei j 



df 



de: 



df 



i/r 



(23) 



In order to take into the account the multiplicity, which results of this change of variables, 
we will divide the range of integration in the e space into domains mi, I = l,...,M , such 
that in each domain the sign of the Jacobian does not change. This results in 

M „ |A r j+ 



(2 A) l=1 Jm H=li¥zj Jfi 



' l , , 1 

Ji df ■ 



i/r 



(24) 



where Ji is the Jacobian 



df 



in the I — th domain and /, are the extremal values of / in 



this domain. Such splitting into a finite number of domains is possible due to Lemma [5] and 



since 



df 



1 

v/2 



de 3 de j+1 



Since the integrand is positive the value of the integral can 
be only increased by increasing the integration volume. Therefore we obtain the following 
inequality 



l/l 



k=l 



<max|^J^|c fc | = (2 + A)^|c fe | <Q 



(25) 



where we have used the fact that the eigenvalues are bounded and the sum is finite. The 
integral can be bounded by 



M 



(2A) l=1 om H=x ^j 



< 



df-Fz sup 



df 




(26) 



dei 



with \Ji\ < 



su P{ £l } 



df 



9el 



-1 



the volume transforms as 



(2A) l=1 Jm H=1 ^ 

The last term is an integral over all Sj excluding one, therefore 

(27) 




Using Corollary [3] and Remark [6] the integral can be further bounded by 

Q 1 - 3 



I f < C7 1 - 



J& 2A(l-s) 

with Ds independent of of the size of system and of R while Cs is the smallest Cj. 



(28) 



IV. THE EFFECT OF RENORMALIZATION 

In some applications 0, [ljj the function / = CkE ik should be replaced by /' = 
^ k C/gEl , where E[ are the renormalized energies and /' = f+/3F({e}). In this application, 
to the leading order only one energy denoted by E is renormalized, namely 

E{ = Ei + (3S h0 V 000 

where (3 is a small parameter and V^ 000 = Xh^BoW- We show that to this order Cs is not 
affected by F for (3 <const e~ 2 ^ 7max_7min ' ) ' :r ' 5 ', where 7 m i n and 7 max are the minimal and the 
maximal Lyapunov exponents and \x$\ is defined by C$ = (Cj,j = xg) . 
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The derivative of Vq is given by 



E 4 ^ 



(29) 



Lemma 8 



Pr 



dV 000 



> D e - 2 (Tmin-^)l^l < o < s < 1 



(30) 



where 7 m i n «s i/ie minimal Lyapunov exponent. 

Proof. Using a perturbation expansion in Ej (to the first order), this derivative, 8 ' 



be calculated 



—q— - v*. u) 2^ En _ Eh 



o^, can 
(31) 



" J k^n 

We now use the fractional moment technique to get a probabilistic bound on this derivative. 
Using the bound on the eigenfunctions and the estimate on the average of ^ E l^l 3 / > ( see 
Theorem [7]), we obtain 



9ipE n (i) 



E 

k^n 



*pE k {j) 1pE k (i) 



E„ — Eu 



<su P [\^ En (j)rEi^^n^ 

w \ k+n 
< £) e £) 5 Q — faito*kj-!»nl e -7jnja*l*j— x i\ 



\E n — Eu 



(32) 



Therefore the derivative 8V q £ (where Xq — 0), can be bounded by 



<9U 000 



< 4D e D$e~' yininS \ Xj \'y ^e~ 37minS ^^e~ 7minS ^ J ~ x ^ 

i 

< 4j) E £) (5 e^ 27minS ^ J iy^e~ 27minS ^^ = D £ Dse~ 2rymluS ^ Xj K 



(33) 



Using the definition of e j 



dV 000 



del 



< 



2 s/2 



dV 000 



dV 000 



de 



i+i 



5/ 



Utilizing the Chebychev inequality we complete the proof. 



Remark 9 Since the derivative fl ° . «s bounded from above we can always select (3 in such 



de 



way that its contribution will be smaller than the derivatives of the energies. 
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